On dynamical tunneling and classical resonances 
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This letter establishes a firm relationship between classical nonlinear resonances and the phe- 
nomenon of dynamical tunneling. It is shown that the classical phase space with its hierarchy of 
resonance islands completely characterizes dynamical tunneling. In particular, it is not important 
to invoke criteria such as the size of the islands and presence or absence of avoided crossings for a 
consistent description of dynamical tunneling in near-integrable systems. 



Dynamical tunneling as a concept emerged more than 
two decades ago in the field of chemical physics where 
it results in the transport of vibrational quanta be- 
tween degenerate modes - a process that would be clas- 
sically forbidden. The importance of dynamical tunnel- 
ing in the molecular context can be hardly overstated 
since this phenomenon provides a route to energy flow 
through the molecule in the absence of direct classi- 
cally resonant mechanisms. Early pioneering worki*2*&^& 
mainly by the chemical physics community provided both 
semiclassicali^ and purely quantum perspectives^ on 
dynamical tunneling. Semiclassically the phase space is 
the natural setting whereas the quantum approach in- 
vokes high order perturbation theory involving a chain of 
off- resonant virtual states (vibrational superexchange 4 ). 
Although seemingly different, there are hints 4 *^ towards 
a connection between the two perspectives and this paper 
attempts to provide further clues. 

The initial suggestion 2 regarding the importance of 
phase space structures to dynamical tunneling has 
been intensely studied and established by the nonlin- 
ear dynamics community over the last decad o 7 ! 8 ! 9 ! 10 ^ 1 ! 12 . 
Dynamical tunneling is found not only to be in- 
fluenced by chaos 8-9,12 but also by various nonlin- 
ear resonances^S* 1 ^* 1 ^ with some recent experimental 
support 13 . In the molecular context Heller recentljii^ 
made a number of interesting observations and conjec- 
tures on the possible implications of dynamical tunneling 
on high resolution molecular spectra™. The most impor- 
tant amongst these is the claim that a nominal 10 _1 -10 -2 
cm -1 broadening of spectroscopically prepared zeroth or- 
der states is due to dynamical tunneling between remote 
regions of phase space facilitated by distant resonances. 
Arguments were provided for identifying the specific res- 
onances and subsequent calculation of the splittings. The 
purpose of this letter is to confirm the above claim via 
a detailed analysis of a relatively simple, albeit realis- 
tic, model spectroscopic Hamiltonian. The analysis also 
indicates that the correspondence between classical res- 
onances and avoided crossings, while interesting, is not 
needed for an understanding of dynamical tunneling. 



We use a model spectroscopic Hamiltonian 1 £ 



H = H + gV}a ] 



■Kv^ + v^) (i) 



with 



H = bJ s {n\ + n 2 ) + u b nb + x s {n\ + n 2 ) 

-I- x b n\ + x sb n b (ni + n 2 ) + x ss n x n 2 (2) 

The values of these parameters, in cm -1 , lu s — 
3885.57, ui b = 1651.72, x s = -81.99, x h = -18.91, x M = 
— 12.17, and x sb = —19.12 are representative of the H2O 
molecule^ 7 ". The three anharmonic modes are labeled as 
stretches (1,2) and a bend (b) and H is symmetric un- 
der 1 <-> 2. The j th mode occupancy is rij = a^aj with 

(aj,Oj) denoting the usual harmonic oscillator creation 
and destruction operators for mode j. The various per- 
turbations Vp:p — (a\) q (a,j) p 4- h.c. connect zeroth-order 



states In), In') with 



q and 



n ; 



The 
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classical limiU^ of the above Hamiltonian is a nonlin- 
ear multiresonant Hamiltonian H(I,9) with (Ij,9j) cor- 
responding to the action-angle variables associated with 
the mode j. In particular the classical limit of V p \q is of 
the form 2^IfI? cos(qdi — p9j). H can be obtained by 

a fit to the high resolution experimental spectra or from 
a perturbative analysis of a high quality ab initio poten- 
tial energy surface. In cither case such effective Hamil- 
tonians provide a very natural and convenient represen- 
tation to understand the spectral patterns of molecular 
systems 18 . Note that despite the three coupled modes, 
H is effectively two dimensional due to the existence of 
the conserved quantity P = (m + n 2 ) + n b /2 called as 
the polyad number. The classical Hamiltonian is inte- 
grable if j3 = and for our choice of parameters it is 
near-integrable if /3 =/= 0. Throughout this study we fix 
P = 8 and (3 — 26.57 cm -1 and denote, due to conserved 
P, the zeroth-order states by \ni,n 2 ). 

To begin with consider the case wherein only the 2:1 
resonances are present i.e., g — 7 = 0. In order to em- 
phasize and illustrate the concept we choose the zeroth- 
order degenerate states |2,0) and |0, 2) without loss of 
generality. Since (2, 0\V^ |0, 2) = = (2,0|t^ 6) |0,2) 
the state |2, 0) is uncoupled from the symmetric counter- 
part |0,2). However, dynamical tunneling can mix these 
states and indeed from Fig. ^ one observes a coherent 
transfer of population with a period of about 0.15 ns 
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FIG. 1: State or zeroth-order number space with the location 
and width of the 2:1 resonance zones. An example superex- 
change path (dashed) between the degenerate states |2, 0) and 
|0, 2} is indicated. The inset shows the survival probability for 
|2,0) and |0, 2) (symbols) versus time. 

corresponding to to a splitting A2 ~ 0.22 cm . The 
nontrivial nature of this process is amplified when one 
considers the fact that at the energy corresponding to 
1 2,0) the primary 2:1 resonances are absent in the clas- 
sical phase space. As shown in Fig. ^ the two states 
are far away in state space from the 2:1 primary reso- 
nance zones. One possible explanation of the tunneling 
arises from the perspective of high order perturbation 
theory or vibrational superexchange^. In this approach 
the states coupled locally by the 2:1 perturbations are 
considered and one constructs perturbative chains which 
connect the two states |2, 0) and |0, 2). An example of 
such a chain is |2,0) -> |1,0) |1,1) -> |0, 1) -> |0,2). 
The contribution to the splitting from the chain is given 
by perturbation theory to be 

4 (20|y( lb )|io)(io| i i/( 2b )|ii)(ii|t>( lb )|oi)(oi|y( 2f ')|02) 

(3) 

with AE° i) „ 2i „ 6 ee E° 012 - El un ^ nb . In principle there 
are an infinite number of chains that connect the two de- 
generate states. In practice, due to the energy denomina- 
tors and near-integrability, it is sufficient to consider the 
minimal length chains^. In our case there are six minimal 
chains and summing the contributions from each one of 
them one obtains a splitting of about 0.23 cm -1 . This 
compares well with the exact splitting but it is impor- 
tant to note that all six perturbative terms have to be 
considered for this agreement. Note that although for- 
mally the superexchange approach invokes the resonant 
terms the connection to the classical phase space is lost. 
If indeed dynamical tunneling is properly understood in 



FIG. 2: Surface of section at E = £2,0,12 and (3 = 26.57 
cm -1 for varying primary 1:1 strength in cm -1 , (a) g — 0, 
(b) g = -2.7, (c) g = -2.95 « go, and (d) g = -3.4. Note 
that for g < go the island corresponds to the induced 1:1 
whereas for g > go the island correspond to the primary 1:1. 
Also note that a Husimi representation of the states |2, 0) and 
|0, 2) would be localized about (I\ — i2)/2 = ±1 respectively. 

the phase space then surely there must be a phase space 
analog of the superexchange approach. The rest of the 
paper is dedicated to uncovering precisely such a phase 
space picture. 

As mentioned the primary 2:1 resonances do not ap- 
pear in the phase space at E — E® 12 and hence a direct 
involvement is ruled out. Nevertheless a weak overlap be- 
tween the 2:1s can result in an induced 1:1 resonance 
which can then mediate dynamical tunneling between 
the states. In Fig. [2^ we show the surface of section at 
E = El 

12 an d one indeed observes a resonance island 
between the two states. In order to confirm the nature of 
this resonance zone we use standard methods of nonlin- 
ear dynamics^ to extract the necessary information. In 
essence one starts with the classical Hamiltonian involv- 
ing only the 2:1 perturbations in the formiS 

H(J,tl>;N) = H (J;N) + ep c (N-2J 1 -2J 2 ) (4) 

X [y / JiCOsV'i + y/j^COS^} 

with p c — /3/V2 and N being the classical analog of 
the polyad number. A formal parameter e has been in- 
troduced with the aim of perturbatively removing the 
2:1 resonances, characterized by "01,2, to 0(e). This 
can be done by invoking the generating function F = 
J\i\)\ + J21P2 + e[gi sin -0i -(- g 2 sin -02] where the functions 
.91,2 = 9i,2{J\, J 2) are determined by the condition of the 
removal of the primary 2:1s to 0(e). The angles con- 
jugate to J are denoted by -0. The procedure is alge- 
braically tedious and we skip the details to provide the 
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important results. The choice of the functions g\ t 2 turns 
out to be: 



91,2(^1, h) = -(3 C 



(N - 2Ji 



fl s + 2a ss J 12 + OLviJiA 



(5) 



where Q, s = lu s — 2u>b + (x s b — Axb)N, a ss = x s + Axb — 2x s b 
and «i2 = Sxb — 4x s b + x ss . Using the above result it is 
possible to show that an induced 1:1 resonance appears 
at 0(e 2 ) with a coefficient 



Qirid 



2 



2J 1 -2J i )f(J 1 ,J 2 ;N)VJxJa (6) 



with f(J\,J 2 ]N) being a complicated function of the 
actions. At this stage the transformed Hamiltonian to 
0(e 2 ) still depends on both the angles "ipx and ij>2 an d 
hence non-integrable. In order to isolate the induced 
1:1 resonance we perform a canonical transformation 
to the variables (K, 0) using the generating function 
G = (4>i — ^ 2 )Ki/2 + + ^2)^2/2 and average the 
resulting Hamiltonian over the fast angle <p2- The reso- 
nance center, K \ — 0, approximation is invoked resulting 
in a pendulum Hamiltonian describing the induced 1:1 
resonance island structure seen in the surface of section 
shown in Fig-E^t. Within the averaged approximation the 
action K2 = J\ + J2 is a constant of the motion and can 
be identified as the 1:1 polyad associated with the sec- 
ondary resonance. The resulting integrable Hamiltonian 
is given by 



H(K 1 , ( p 1 ;K 2 ,N) 
where 

9ind(K 2 ,N) 

with 

f(K 2 ,N) = 



1 



2Mi 



-Kl + 2g ind (K 2 ,N) cos 20i 



2(aia - 2a ss )~ 1 



^-f(K 2 ,N)(N-2K 2 )K 2 



A(n s +a ss K 2 ) + ai 2 N 
[2(O s + a ss K 2 ) +a l2 K 2 f 



(J) 



(8b) 



(8c) 



In terms of the zeroth-order quantum numbers K\ = 
n\ — n 2 and K2 = n\ + ri2 + 1 = m + 1. 

One can now use the above pendulum Hamiltonian to 
calculate the resulting dynamical tunnel splitting of the 
degenerate modes \n\ = r, 112 = 0,rib = 2{P — r)) and 
\rii = 0, n 2 = r,n b = 2{P - r)) via4 



Qind 



(r-2) 

n 

=-(r-2) 



9i. 



E* R {r) - E° R (m) 



(9) 



where E R (k) = k 2 /2Mn is the zeroth-order energy. For 
our example with r = 2, to = 2 using the parameters 
of the Hamiltonian we find Mn w 1.32 x 10~ 2 and 




4.43 cm 1 . The resulting splitting Aj 



0.26 



FIG. 3: The variation of the dynamical splitting A2 with the 
primary 1:1 resonance strength g is shown (solid line). The 
WKB estimate is shown as dashed line. WKB estimate with 
a small 2:2 resonance added is also shown (circles). The inset 
shows a similar calculation with a different set of parameters 
representing the D2O molecule—. In this case (3 = 15.78 and, 
g inc i ~ 2.73 cm -1 . 



cm -1 agrees very well with the exact splitting. This 
proves that the induced 1:1 resonance arising from the 
interaction of the two primary 2:1 resonances is mediat- 
ing dynamical tunneling between the degenerate states. 
At this juncture it is important to note that the induced 
resonance strength is quite small and the two states are 
not involved in any avoided crossing. Moreover, from a 
superexchange perspective it is illuminating to note that 
the splitting can be calculated trivially by recognizing the 
secondary phase space structure (a viewpoint emphasized 
in Ref. [ll| as well). In comparison the original superex- 
change calculation, without any reference to the phase 
space, required taking into account 6 terms with varying 
signsSS. This observation emphasizes the superior nature 
of a phase space viewpoint on dynamical tunneling. 

As a further demonstration of the role of nonlinear res- 
onances in dynamical tunneling we consider the Hamil- 
tonian in eq. ^ with 5^0 with 7 = 0. In particular the 
sign of the primary 1:1 perturbation strength g is taken 
to be the opposite of the induced 1:1 strength j m( j. If the 
induced resonance is playing a role then from our anal- 
ysis we expect that the primary and induced resonances 
will cancel each other around g = go ~ —2gi n d/(m + 1) 
resulting in small splittings in this region. In Fig. |3J the 
exact splittings are shown as a function of g with the 
WKB results^ for comparison. This confirms our ex- 
pectations to a certain degree in that the splittings are 
undergoing dramatic changes in the vicinity of go. A 
crucial observation is that the exact splitting is orders 
of magnitude larger than the simple WKB estimate and 
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become small slightly away from go- On the other hand 
the classical phase space in Fig. ^indicates the predicted 
disappearance of the 1:1 islands. From our arguments 
this far it would be natural to associate one or more high 
order nonlinear resonances with the residual tunneling 
around go since the simple semiclassical estimate for go 
was based on the 0(e 2 ) induced 1:1 resonance cancelling 
the primary 1:1 resonance. In reality there are the har- 
monics of the 1:1 resonance that appear at higher orders 
in e. It is expected that the strengths of such higher 
harmonics like 2:2, 3:3, etc. would be extremely small. 
Nevertheless around g — go the most dominant resonance 
involved in dynamical tunneling would be the 2:2. The 
strength of this tiny but dominant 2:2 resonance can be 
estimated roughly by adding a 2:2 perturbation (7) to 
eq. ^with (3 — 26.57 cm -1 , g = go and noting the value 
of 7 for which the exact and WKB results come close. 
A much more rigorous estimate, which is a difficult cx- 
cersise in classical perturbation theory, can be made by 
going to higher orders, atleast 0(e 3 ), in e. We now esti- 
mate the splitting with a WKB calculation including the 
2:2 resonance with strength 7 w 2.32355 x 10~ 4 cm -1 . It 
is clear from Fig. [3] that the exact splitting and the mod- 
ified WKB estimate based on the higher order 2:2 agree 
fairly well. It is also satisfying to see that the modified 
WKB calculation hardly effects the splittings far away 
from g — go- As an independent check in Fig. [21 inset) 
we show the same calculation for a different set of Ho pa- 
rameters representative of the D2O molecule^ and the 
results are similar. This supports the argument that in 
the vicinity of go, where the 1:1 resonance is absent, the 
extremely small 2:2 resonance is mediating the dynam- 
ical tunneling between the states |2,0) and |0,2). Two 
remarks are in order at this stage. First the two states do 
not undergo any avoided crossing as a function of the pa- 
rameter g. This can also be indirectly inferred from the 
fact that a superexchange calculation of the splittings 
essentially reproduces the exact result and thus include 



the importance of higher order resonances near g = go- 
A more detailed analysis of the perturbative chains from 
the semiclassical viewpoint would be interesting. Sec- 
ond the modified WKB calculation is in good agreement 
with the exact splittings only in the vicinity of go by ne- 
cessity. There are contributions from even higher order 
resonances which are absent from our simplified analysis 
and a subtle interplay of all the nonlinear resonances give 
rise to the exact result. 

To summarize, in this work using a model spectro- 
scopic Hamiltonian we have demonstrated the intimate 
connection between dynamical tunneling and the reso- 
nance structure of the classical phase space. Thus dy- 
namical tunneling connects two degenerate states as long 
as there is a nonlinear resonance juxtaposed between 
them as viewed in the phase space. The order and width 
of the resonance are immaterial. This supports an earlier 
claim regarding the possibility of dynamical tunneling 
as a source of narrow spectral clusters associated with 
spectroscopically prepared, localized, zero-order states. 
However the notion that such resonances are the cause 
of avoided crossings does not seem to hold. Consequently 
it is also not necessary that only a specific classical reso- 
nance be the agent of dynamical tunneling. Primary, in- 
duced and even higher harmonics of the resonances can 
mediate dynamical tunneling and the consequences for 
energy flow and control from this standpoint seem crucial 
and needs further study. It is interesting to note that in 
multidimensional near-integrable systems nonlinear res- 
onances would be involved in two long time phenomena - 
dynamical tunneling and Arnol'd diffusioniS. The com- 
petition between them and their spectral consequences 
are worth investigating from a fundamental standpoint 21 . 
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